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ABSTRACT 

More than 50 years after the dawn of X-ray astronomy, the dynamical parameters of the 
prototypical X-ray binary Sco X-l are still unknown. We combine a Monte Carlo analysis, 
which includes all the previously known orbital parameters of the system, along with the K- 
correction to set dynamical constraints to the masses of the compact object (Mi < 1.73 Mg) 
and the companion star (0.28 M 0 < M 2 < 0.70 M 0 ). For the case of a canonical neutron 
star mass of Mi ~ 1.4 M 0 , the orbital inclination is found to be lower than 40°. We also 
present the best near-infrared spectrum of the source to date. There is no evidence of donor 
star features on it, but we are able to constrain the veiling factor as a function of the spectral 
type of the secondary star. The combination of both techniques restricts the spectral type of 
the donor to be later than K4 and luminosity class IV. It also constrains the contribution of the 
companion light to the infrared emission of Sco X-l to be lower than 33%. This implies that 
the accretion related luminosity of the system in the K band is larger than ~ 4 x 10 35 erg s _1 . 
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1 INTRODUCTION 

Low-mass X-ray binaries (LMXBs) harbour a low-mass donor star 
which transfers matter onto a black hole or a neutron star (NS) via 
an accretion disc. They can be divided into two populations accord¬ 
ing to their long-term behaviour. Transient systems spend most part 
of their lives in a faint, quiescent state, but show occasional out¬ 
burst, where their X-ray luminosity increases above ~ 10% of the 
Eddington luminosity (I/Edd)- There is also a population of ~ 200 
sources that are always X-ra y active , displaying l umi nositie s above 
Lx — 10 36 ergs _1 (but see I Armas Padilla, Degenaar & Wiinanda 
120131) . The so-called persistent sources mostly harbour NS accre- 
tors, whilst the vast majority of black holes are found in transient 
systems. Transient LMXBs have provided a number of dynamical 
NS and black hole mass measurements thanks to th e detecti on of 
the d onor star during the quiescent phase (e.g. ICasares & Jonkeij 
|20I4|) . This is not the case for the majority of the persistent pop¬ 
ulation, where the companion spectrum is totally swamped by the 
reprocessed light from the accretion flow. Only for a few long or¬ 
bital period systems a full orbital solution exists, since they harbour 
giant (i.e. bright), companion stars. 

Due to the blue spectrum of the accretion discs, and the typi¬ 
cally late spectral types of the donor stars, near-infrared (NIR) ob¬ 


servations offer a go od opportunity to detect the companion star 
spectral features (e.g. lBandvopadhyav et alJ 19971) and obtain a full 
dynamical solution (see lSteeghs et alj2013l for astudv purely based 
on NIR data). 

Scorpius X-l is the prototype LMXB and also the bright¬ 
est persistent X-ray source in the sky, being the ta rget o f nu¬ 
merous studies since its discovery IGiacconi et alJ 1 1 9621) . Al¬ 
though the spectral classification of the companion star is un- 
kn own, its r elatively long o rbital period (P OT b = 18.9 h, 

GpttliebjJVrigh^&^ifle^l^T^) sugg ests a n evolved, late type star. 
Br adshaw . Fo malo n t & Geldzahler j 19991 ) measured the trigono¬ 
metric parallax of Sco X-l using VLBA radio observations, and 
deduced a distance of 2.8 ± 0.3 kpc. Further observations allowed 
the detection of twin radio lobes, leading to an inclination value 
of 44° ± 6°. This value assumes that t he rad io-je t is perpendicu- 
lar to the orbital plane dFomalont. Geldzahler & Bradshawll200lh . 
Even though the system does not show any of the NS identifying 
features (thermonuclear X-ray bursts or pulsations), a N S accreto r 
is widely assumed based on its X -ray behaviour (see Ivan der Klisl 
1 20061 : iMunoz-Darias et al1l2014l for X-ray emission from NS in 
LMXBs)._ 

ISteeghs & Casares! d2002l) presented an optical spectroscopic 
technique for measuring system parameters in persistent LMXBs. 
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2 D. Mata Sanchez et al. 


It is based on the discovery of narrow emission lines within the 
Bowen blend, arising from the irradiated face of the donor star 
in Sco X-l and powered by fluorescence. This claim is based 
on the narrowness, Doppler velocities and phase of the lines. 
Also, timming studies of the fluorescence emission have shown 
a time lag with the irradiating X-ray flux which is consistent 
with the light time betwe en the two components of the binary 
(Munoz-D arias et al .1120071) . Radial velocity studies of the Bowen 
lines (Bowen technique here inafter) have been successfully applied 
to oth er dozen LMXBs (see lComelisse et ah 2008; Munoz-Dariasl 
1 20091) . providing the first dynamical solutions to some clas sical 
NS systems (e.g. ICasares et al.ll2006l : IComelisse et al ■112 0071) and 
the canonical black hole transient GX 339-4 l lHvnes et al.ll2003l : 
iMunoz-Darias et al]|2008h . The velocity inferred from these emis¬ 
sion lines (K e m ) corresponds to that of the irradiated region of 
the companion; hence it only represents a lower limit to that of 
its centre of mass ( K 2 )• In order to correct K eul from this effect 
numerical solutions _should be applied, such as those computed by 
Mun oz-D arias . Ca sares & Martfne z-Pai s J2005I . see Sec.Q. 

In this paper we aim at constraining the masses of the com¬ 
ponents in Sco X-l by combining two different techniques: (i) a 
K-correction along with Monte Carlo analysis accounting for con¬ 
strains to the orbital parameters from previous studies (Section [3 
and (ii) NIR intermediate-resolution spectroscopy obtained around 
orbital phase zero, when the non-heated face of the donor star is 
oriented towards the Earth (Section^- We note that Sco X-l, be¬ 
sides being the prototypical LMXB, is also an important object in 
the search for persistent gravitational waves, which al so requires a n 
accurate determination of i ts system parameters (e.g. lAbadie et al .1 
l201ll;lGallowav et al .112014) . 



q 


Figure 1. K-correction for Sco X-l obtained through a Monte Carlo simula¬ 
tion (10® events). We plot in solid and dashed, coloured lines the obtained 
limits to Mi as a function of q. The upper limits correspond to K c (ao) 
(blue, dashed line) and RT max (green, dashed line). The lower limits are 
defined by K c {oim ) (red, solid line) and K\ m i n (black, solid line). The 
yellow band represents the 90% probability region delimited by the previ¬ 
ous limits. The violet, dotted line corresponds to the imposed Mi > 1 M@ 
constraint (see text), which restrict the allowed parameter space to the or¬ 
ange, stripped triangle. 


0°). This gives the maximum K-correction and, therefore, K c 
would be minimum. 


2 THE K-CORRECTION 


The companion star of Sco X-l has been solely detected through 
the Bowen t echnique described a bove. A K em velocity was initially 
reported by Isteeg hs & Casaresl l2002i ). and very recently refined 
bv lGallowav et alj 1 2014 , see Table QJ. This velocity can be trans¬ 
lated to that of t he centre of mass of the don or by the so-called K- 
correction (K c ). IMunoz-Darias et all d2005l) presented numerical 
solutions for K c , which is parametrized by the following equation: 


K c = 


Ke , 


K 2 


= N 0 + N 3 q + N 2 q + N 3 q 3 + A r 4 q 4 


0 ) 


where q = is the mass ratio (Mi and M 2 are the masses of 
the compact object and the donor star, respectively). The Ni values 
are tabulated and depend on both the opening angle of the accretion 
disc that shades the companion (a), and the orbital inclination ( i ). 
Therefore, K c is always lower than 1, and approaches unity as a 
increases and the regions of the companion with lower radial ve¬ 
locities become shaded. For the limit case in which the companion 
is fully shaded by the disc (i.e. a = <xm), K c is at maximum and 
can be approximated by the following expression. 


K c * 1 - f(q) • (1 + q) 


/(«) = 1-0.213 



( 2 ) 


Therefore, for a given q and K e m , the true K 2 velocity is in 
the range: 


K eul /K c (ctM) < K 2 < K eul /K c (ao) 


(3) 


where K c (a 0 ) represents the case in which the disc is not present 
and the inner face of the companion star is fully irradiated (a = 


2.1 Monte Carlo analysis 

The mass of the compact object (Mi) can be inferred using the 
mass function formula, which is directly derived from Kepler’s 
laws: 


= fU Al(l + q ) 2 
2-kG sin 3 i 


(4) 


For the case of Sco X-l, some of the above system parameters are 
poorly constrained (e.g. orbital inclination) and for some others we 
just have extremal values (see also Section |4]). This is the case of 
AT and K 2 , the radial velocities of the compact object and the 
donor star, respectively. To deal with this, we performed a Monte 
Carlo analysis in order to obtain reliable limits to the compact ob¬ 
ject mass. We simulated 10® random, normally-distributed values 
of i and K e m , and lower and upper limits to the radial velocity of 
the compact object (K i m in and ATmax, respectively) centred at the 
values reported in TableQ] We apply both K c (ao) and A'c(qm) to 
each simulated K em value. Finally, we solved Equation [4] for 10 2 
steps in q (from 0 to 0.8), using the previously defined distribution 
o f i va lues, as well as the orbital period reported in lGallowav et al] 
J2014t) . We obtained four M 1 distributions: two lower limits cor¬ 
responding to K 2 (aM) and ATmin, and two upper limits from 
K 2 (aq) and ATmax. These limits are represented by solid lines 
in Fig. Q] and are obtained by adopting the value of Mi above or 
below 90% of the points (for upper and lower limits, respectively) 
in each distribution and for every sampled value of q. 
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Figure 2. Normalized NIR spectrum of Sco X-l obtained by averaging the 7 nodding cycles (black soli d line). For comparison we s how (blue solid line; 
positive offset applied) the spectrum degraded to the resolution of the best NIR spectrum available so far iBandyopadhyay et alJI 1997b . We also include the 
spectra of three template stars (negative offsets applied) from IRTF Spectral Library as black and green solid lines. The latter have been veiled by the limit 
values reported in Table [2]for its corresponding spectral type. Typical NIR spectral features are identified by vertical, red, dashed lines. 


Table 1. Observational measurements included in the Monte Carlo analysis. 


ftrn[kms x ] .fs'lmiljkms J ] ATmax[kms *] i[deg] 


1201 3h . we further restrict the mass ratio to 0.28 < q < 0.51, and 
hence 0.28 Mg < M 2 < 0.70 Mg. 


74.9±0.5« 40 ±5< 2 > 53 ± l( 2 > 44 ± 6® 3 NEAR INFRARED SPECTROSCOPY 


R EFERENCES: (l)lGallowav et all 120141) : t2l lSteeghs & Casared J2002h : 
I'f l ll omalont et al.l 1200 lh : 


2.2 The masses of the NS and the companion 

The intersections of the solid lines in Fig.Q]provide absolute limits 
to Mi and q. They yield extremal values to the system parameters 
with, at leasfj 90% confidence level: 


0.22 Mq < Mi < 1.73 M 0 
0.24 < q < 0.75 
0.05 M 0 < M 2 < 1.30 M 0 


therefore, we obtain a robust upper limit of Mi < 1.73 Mq, con¬ 
sistent with a NS accretor. The lower limit is, on the other hand, 
not very restrictive. Adopting instead Mi m i n = 1 M q as a conser- 
vative minimum value for the mass of a NS (see e.g. Kiziltan et al.1 


We observed Sco X-l on 2008 June 28, using the Long- 
slit Intermediate Resolution Infrared Spect rograph (LIRIS, 
Acosta-Pulido et alj2003l .rMan chado et al.l2004l) at the Cassegrain 
focus of the 4.2-m William Herschel Telescope (WHT), in La 
Palma (Spain). We used the K-band grism (R=2500, 3.5 A/pix dis¬ 
persion) combined with a 0.75” slit. The use of grisms as dispersers 
usually does not permit to vary the spectral range. In our case, the 
employed slit (10p75ext) consists of 2 slits which are offset along 
the dispersion axis by about 115 pixels in both senses. This allowed 
us to redshift the nominal spectral range to 2.09 — 2.45 gm, which 
includes the CO bandheads. 

We performed seven nodding cycles of 3 positions each, with 
an exposure time of 900 s per cycle at airmasses in the range 
~ 1.4 — 1.5. We cover orbital phase s, if = 0.01 — 0.13 (using 
the ephemerides from lGallowav et ahll2014l) and the total exposure 
time was 1.75 h. For the data reduction process we used IRAf0 


1 The confidence level of these results is actually slightly higher since two 
90 % limits (apart from that of the inclination) are simultaneously verified. 


2 IRAF is distributed by National Optical Astronomy Observatories, op¬ 
erated by the Association of Universities for Research in Astronomy, Inc., 
under contract with the National Science Foundation. 
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lirisdr task packag^H as well as a modified vers ion of the package 
xtellcor_general developed by lVacca et all20 03l to perform the tel- 
luric correction to the final spectra (see lRamos* Almeida et al.l2009l . 
120131 for a more detailed description of the data reduction process). 

This strategy yielded seven (nodding averaged) NIR spectra 
of Sco X-l. Absorption features from the donor star are not ev¬ 
ident in any of them, which were combined to produce a higher 
signal-to-noise average spectrum. Points deviating more than 3 a 
above or below the continuum, probably caused by bad pixels or 
residual from the sky subtraction, were interpolated. The com¬ 
bined, normalized spectrum is presented in Fig. [2] To our knowl- 
edge this is the best NIR spe ctrum of Sco X-l available to date (see 
iBandvopadhvav et al.l 1 19971 for a previous example). No absorp¬ 
tion lines from the companion star are detected. The most promi¬ 
nent feature in the spectrum is a strong Br 7 emission line, whose 
full-width-at-half-maximum, FWHM = 480 ± 30km/s, is con¬ 
sistent with that of Ha (FWHM = 440 ± 20km/s). This latter 
value was obtained from previous observations performed by our 
team with the Intermediate Dispersion Spectrograph (IDS) at the 
Cassegrain focus of the 2.5-m Isaac Newton Telescope (INT), in 
La Palma (Spain). 


3.1 Veiling factor 


Since no features from the companion are evident in our NIR spec¬ 
trum, we investigated the amount of flux from the accretion flow 
required to veil the donor. This obviously depends on its spectral 
type. As a first step, we compared our spectrum with templates 
of G-M main sequence stars from the IRTF Spectral Librar 
In Fig. [2] we present normalized spectra of GOV, KOV and M5V 
stars, which cover the same spectral region as that of Sco X-l. For 
each spectral type, we define Fd ep th as the normalized flux of the 
deepest photospheric absorption line present in the NIR spectrum. 
These correspond to Nal (2.206 pm), Mgl (2.282 pm) and 12 CO 
(2.295 pm), for M5V, GOV and KOV, respectively. Similarly, we 
define the veiling factor ( X ) as the fractional contribution of the 
accretion related luminosity (L^ cc ) to the total flux. The veiling 
necessary to make these features shallower than the noise level (3<r) 
within the corresponding spectral region of the Sco X-l spectrum 
is a lower limit to this factor. Therefore, 


V > 1 


3(7 

Fdepth 


X = 


t k 

- LJ acc 

I I< 1 tK 

J-'acc 1 J-s 2 


(5) 


where Lf' is the luminosity of the donor. X is always in the range 
0 X < 1 and depends on the spectral type of the companion. 
Our lower limits to the fractional contribution of the accretion lumi¬ 
nosity of Sco X-l NIR flux (A' m in) are reported in the first column 
of Table [2] It varies from 0.09 to 0.78 for the cases of GOV and 
M5V, respectively. 


4 DISCUSSION 

The K-correction and Monte Carlo analysis (see Sec. [2} im¬ 
pose an upper limit to the compact object mass of Mi < 
1.73 Mq. This rules out a massive neutron star (~ 2 Mq) like 
e.g. that presented in [Pemorest et al .1 j201 (I) . More accurate val¬ 
ues of a, i and A'i are necessary to improve this constraint. 


3 The IRAF lirisdr task package is supported by Jose Acosta. 

4 http://irtfweb.ifa.hawaii.edu/~spex 


Ide Jong, van Paradiis & Augusteiinl i l 19961) proposed an average 
value of a = 12° for persistent LMXBs. Using this value, we can 
calculate K 2 by applying the same methodology, which constrains 
the mass of the compact object to Mi < 1.25 Mq. This value is 
lower than the canonical NS mass (Mi ~ 1.4 Mq). However, we 
note th at NS masses d own to Mi ~ 1.2 Mq have been reported 
(e.g. iKiziltan et al]l2013l) . If we fix the NS mass to the canonical 
value, we obtain a ^ 11° and K\ > 47 km s -1 . All the above 
constraints require an inclination as low as i ~ 36° (we are quot¬ 
ing 90% confidence limits). 

The upper limit to K\ used in this work was obtained by 
applying the double Gaussian method to the Hell lines, while 
the lower limit comes fro m Doppler tomography analysis of H/3 
dSteeghs & Casaresll2002l) . Our analysis rules out this lower limit 
for a NS mass higher than 1 Mq (see FigO). and suggest the upper 
limit to be close to the real value. Interestingly, the latter stands true 
for another NS system analysed using the double gaussian tech¬ 


nique, namely XI822-371 (Munoz-Darias et al. in prep.; see also 
|Casares_gt_alJ 2003 ), in whic h the Ki velocity is accurately known 
(J onker & van der KlisboQll) . Fixing Mi = 1.4 Mq and assuming 
K 1 = Aimax, we obtain an absolute upper limit of i ~ 40°, which 
is only satisfied if a = 0° (i.e. applying K c (a 0 )). Larger values of 
the opening angle or lower values of K\ would result in even lower 
inclinations. We note that inclinations > 50° were favoured to bet¬ 
ter explain time lags between reprocessed light (Bowen emission) 
and X-rays observed in Sco X-1. This value was found when the 
time-lags were compared with those predicted by numerical trans¬ 
fer functions with a maximum response. However, we note results 
consistent with those presented here are also obtained if a broader 
range of expec ted delays and/or more co nservative solutions are 
considered (see lMunoz-Darias et al.ll2007l for details). 

Finally, since Sco X-l is not eclipsing, we can use our con¬ 
straints to the mass ratio (q) and the Roche lobe radius (derived 
through the formulae in IPaczvriskilll97lh to infer a maximum in¬ 
clination purely based on this condition (i.e. absence of eclipses). 
This yields i < 70°. 


4.1 The nature of the donor star 

Although our observations were performed at the most favourable 
orbital phase, we fail to detect any companion star feature in the 
NIR spectrum of Sco X-l (see Sec. [3}. Nevertheless, we obtain a 
lower limit to the veiling factor as a function of the spectral type 
of the secondary, and this can be used to discuss its nature. We 
initially found 0.05 Mq < M 2 < 1.30 Mq, and assuming an ac- 
cretor h eavie r than 1 Mq, we infer 0.28 Mq < M 2 < 0.70 Mq. 
Faulkner, Flannery & Warner! 119721 ) presented a direct relation be¬ 
tween the mean density of the companion star and the orbital period 
in the Roche lobe filling binaries. For a main sequence companion 
in a 18.9 h orbital period, this result in a spectral type A0 l iCoxl 
120001 ) . However, the features of an A0V donor would dominate the 
observed spectrum. Similarly, its mass (M2 = 2.9 Mq) is ruled 
out by our results. This implies that the donor star is not a regular 
dwarf star but a somewhat more evolved object. Our upper limit 
(M 2 < 0.70Mq) contraints the spectral type to be later than K4, 
since this value corresponds to a main sequence donor. Using the 
same procedure as in Section [3] for a K4 donor we obtain a lower 
limit to the veiling factor of X > 0.66. 

We can independently test the allowed spectral types for the 
companion star by solely using the veiling factor constraints. In 
a first step, we use the abso lute K -magn itude value of Sco X-l 
M“ = —1.14 (see IWachter et al.ll2005l and references therein) 
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Table 2. Minimum veiling factor (X m i n ) and implied minimum K- 
magnitude (Mk m i n ) for the donor star. M 1 ^ and M^ 1 refer K- 
magnitudes of dwarf-to-giant standard stars. 


Sp. Type 


% min 

MV 

M W 

MW 

GO 

0.09 

-1.04 

a.isC 1 - 2 ) 

3.O5U’ 2 ’ 3 ) 

—1.18( 1,2 ) 

K0 

0.61 

-0.12 

4.13( 1 ’ 2 > 

3.89(T 2 . 3 ) 

—1.65( 1,2 ) 

M5 

0.78 

0.51 

7.80( 1 ’ 2 ) 

7.06( 1,2,3 ) 

—6.50( 1,2 ) 


REFE RENCES: (ll lCoiilfeOOOh : (2) lKoomneeil<1983l) : Ol lPaczvriskil 
<197 ll) 


to set a minimum value of the absolute K-magnitude of the donor 
star as a function of the spectral type (see TableO. Then , we c om- 
bine the values tabulated for Mv and (V — K) Isee lCoxIbOQOl and 
iKoornneetll 1 9831 . respectively) to infer absolute K-magnitudes for 
different spectral types and luminosity classes. We conclude that 
all the inspected spectral types would be compatible with the veil¬ 
ing constraints if the companion were in the main sequence ( Mk). 
However, doing the same for giant stars ( M™) we are able to rule 
out luminosity class III (and earlier) for the donor. On the other 
hand, assu ming that the companion star radius is at the Roche lobe 
( IPaczvhskill 19711 ). we obtain subgiant K-magnitudes (MW ) com¬ 
patible with our veiling constraints (see Tabled. 


We also presented the best NIR spectrum of the source to date. 
Despite the non detection of donor star features, our deep observa¬ 
tions constrain the spectral type of the donor to be later than K4 
and luminosity class IV. We obtain a lower limit to the veiling fac¬ 
tor in the NIR of X > 0.66, implying that the accretion related 
luminosity in the K band is larger than a few times 10 35 erg s” 1 . 
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4.2 Accretion luminosity in the NIR 

The above constraint to the veiling factor (X > 0.66) implies that 
the accretion luminosity (L^ cc ) accounts for more than 2/3 of the 
observed K band flux. If we c ons ider the K band luminosity of 
Sco X-1 reported bv lWachter et al.l ( l2005h . and propagate the errors 
on the interstellar extinction and the distance, we obtain 8 ± 2 x 
10 35 erg s -1 . As noted by these authors, the error is not dominated 
by uncertainties in the measurements but reflects real variability 
(~ 25%) in the K band. According to the MAXI all sky monitor 
i Matsuoka et alj|2009l) . the day of our observations Sco X-l flux 
was within 10% of its long-term average. Since data were obtained 
when the donor contribution is at maximum, we can use the limit 
X > 0.66 to constrain the accretion luminosity of Sco X-l in the 
NIR, which is therefore in range 4 x 10 35 erg s _1 < L^ cc < 10 36 
erg s _1 . 


5 CONCLUSIONS 

We performed a Monte Carlo analysis along with the K-correction 
to constrain the dynamical parameters of the prototypical X-ray bi¬ 
nary Scorpious X-l. We obtain the following constraints, at 90% 
confidence level: 

Mi < 1.73 M e 
0.28 < q < 0.51 
0.28 Mq < M 2 < 0.70 Mg 

Thus, the presence of a massive neutron star (~ 2 Mq) in this sys¬ 
tem is ruled out. If we consider standard values for LMXBs, a pos¬ 
sible set of parameters would be a canonical NS of Mi ~ 1.4 Mq, 
an orbital inclination of i ~ 36° and a disc opening angle of ~ 11°. 
Assuming Mi = 1.4 Mq necessarily implies i <> 40°. Higher val¬ 
ues of the inclination would only be possible if the neutron star is 
lighter than the canonical value or any of the measurements used in 
this work is not correct. 
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